The interplay of microscopic and mesoscopic structure in complex networks 



Jorg Reichardt/'^ Roberto Alamino,^ and David Saad^ 

^Institute for Theoretical Physics, University of WOrzburg, WiJrzburg, Germany 
^Complexity Sciences Center, UC Davis, Davis, USA 
^The Nonlinearity and Complexity Research Group, Aston University, Birmingham, UK 

Not all nodes in a network are created equal. Differences and similarities exist at both individual node 
and group levels. Disentangling single node from group properties is crucial for network modeling and 
structural inference. Based on unbiased generative probabilistic exponential random graph models and 
employing distributive message passing techniques, we present an efficient algorithm that allows one to 
separate the contributions of individual nodes and groups of nodes to the network structure. This leads 
to improved detection accuracy of latent class structure in real world data sets compared to models that 
focus on group structure alone. Furthermore, the inclusion of hitherto neglected group specific effects 
in models used to assess the statistical significance of small subgraph (motif) distributions in networks 
may be sufficient to explain most of the observed statistics. We show the predictive power of such 
generative models in forecasting putative gene-disease associations in the Online Mendelian Inheritance 
in Man (OMIM) database. The approach is suitable for both directed and undirected uni-partite as well 
as for bipartite networks. 

Networks are fascinating objects. Charting the interactions between system constituents, abstracted as edges and nodes, has 
allowed us to marvel the interconnectedness of systems and appreciate their complexity. Whether in understanding foodwebs [IJ, 
social communities [2], protein-interaction [3 4j, metabolism [5], neural networks [6J or communication [7] the network- metaphor 
has been highly successful in advancing our understanding of complex systems. While being esthetic charts of complex systems, 
networks only reveal insights through rigorous statistical analysis and modeling. The abstraction of complex systems as networks, i.e. 
graphs connecting nodes through links, provides a tremendous simplification. Not all network constituents are created equal; network 
structure strongly depends on node properties and functions which may differ radically among individual nodes, or may be shared by a 
whole group of nodes. However, in many circumstances, node properties or functions may not be accurately or not completely known 
in contrast to the interactions between nodes, which often can be mapped accurately and efficiently. A primary goal of network 
research is to deduce unobserved, or latent, node properties through structural analysis. 

Methodologically, structural analysis can proceed using either purely descriptive statistics or a generative model coupled with 
inference techniques; both of which are complementary. However, a descriptive approach can only make assertions about a single 
network, while generative models gives rise to a whole ensemble of networks that are statistically equivalent to the observed dataset. 
Such ensembles are vital to the study of dynamical processes on networks but above all, they allow us to potentially differentiate 
between more and less important structural features. We will hence use generative models in our approach. 

Thematically, research on network structural features has two foci. One is the study of microscopic properties, such as very small 
subgraphs of a given size known as motifs [8J and their respective distributions in networks, or the properties of single nodes, like node 
degree, centrality, betweenness or page-rank, which are correlated with latent node characteristics such as expansiveness, esteem, 
functional importance or authority. The second investigates mesoscopic structural features indicative of properties or functions 
shared by groups of nodes. Differing link densities between entire classes of nodes, assortative or dissortative mixing [9J, community 
structure [10 llj or more generally, block structures [12| \13\, are all examples of the latter. Interestingly, the two themes run 
in parallel. When modeling degree distributions or analyzing motif distributions, group effects are rarely taken into account, while 
individual node properties are generally neglected in inferring latent node classes from network structure via community or block 
structure detection algorithms. 

With this paper, we intend to close this gap. We present a principled probabilistic approach to the inference of latent node classes 
based on a generative model that includes node specific features; these provide a more realistic model and enable one to quantify 
probabilistic measures and confidence levels of the observed structural details. To estimate model parameters, we employ distributive 
message-passing techniques, with computational complexity scaling linearly with the problem size. Using real world data, social and 
biological networks, we demonstrate the significance of node specific parameters to the quality of latent class inference. Through 
an example from neuroscience, we show that including group specific effects in the random null model used to assess the statistical 
significance of motif counts may provide a group based explanation for many motifs. Finally, we demonstrate how such generative 
probabilistic models can be used for the prediction of putative new gene-disease associations. 



I. METHODS 

A. Exponential random graphs 

When inferring structural features in networks with probabilistic generative models, we have to establish what constitutes a good 
model. Within the infinite set of possible models, exponential random graph models (ERGMs) pJjilSj stand out as combining desired 
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FIG. 1: Factor graphs and an example of an elementary message passing update. Factors of the likelihood function are 
represented as squares, variables of the generative model as circles. Connections indicate which variables enter the calculation 
of which factor. A) For a bipartite actor-event networks represented by an N x M adjacency matrix Ai^, class label ai and 
activity ai of actor i enter in the calculation of all factors in row i. Equivalent ly, class label T/j^ and popularity of event 
fi enter in the calculation of all factors in column fi. The variables Brs denoting preference of actors in class r for events in 
class s enter in every factor. Note that while each factor depends on only 0(1) variables, the a and a variables enter in the 
calculation of 0(N), the r and /3 variables in 0{M) and the Brs variables in 0(NM) factors. B) Pictorial representation of 
the messages involved in calculating Ri^(ai) sent from factor Ai^ to variable ai according to equation C) For directed 
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properties that make them the preferred choice for the task. ERGMs are mean unbiased and make the observed data maximally likely. 
They are maximum entropy models thus ensuring that one avoids adding redundant assumptions to the network's given structural 
features. In other words, they parameterize the largest ensemble of networks compatible with our observations, while making the 
observed network typical for the ensemble. Plus, their parameters have a very intuitive interpretation. 

Within the framework of ERGMs it is easy to combine node specific with group specific effects. Consider a given, bipartite network 
specified by \ts N x M adjacency matrix A, representing for instance the attendance of N actors in M events. If actor i has attended 
event /i, then Ai^ = 1 and otherwise Ai^ = 0. Equally, A could represent the choices of N consumers from a list of M products, or 
the association of diseases with M different genes. The possibilities are many and we will use the actor-event picture, but without 
limiting the applicability of the model to this case alone. 

We restrict ourselves to dyadic models, i.e. we assume the entries of the adjacency matrix to be modeled by the conditionally 
independent random variables Di^ e {0? 1}- A simple ERGM that captures both individual (actor- and event-specific) and group- 
specific factors is given in terms of the odds ratio of actor i attending event p [44]: 

= l|g) _ at Bg.r^ 
V{Di^=0\e) (1-aO (1-/3^) (1-B..rj' ^' 

The shorthand denotes the set of all model parameters, i.e. in this case = (ai, .., aAr, .., /3m, cri, .., cta^, n, .., tm, B). Of these, 
only a small subset is relevant for an individual dyad Di^. The parameter ai e (0,1) denotes the global activity of actor i, higher 
ai meaning higher odds of attending any event. Correspondingly, e (0, 1) denotes the global popularity of event p. Furthermore, 
every actor i and every event p carry a class index ai and r^^x, respectively [45]. The matrix Brs ^ (0,1)V r, s, models the data at a 
coarser, group specific level, denoting the tendency or preference of an actor of class r to attend an event of class s. Higher entries 
mean higher odds for the attendance of any actor of class r to any event of class s. The matrix Brs is also called a block model of 
the data. 

The rich literature on ERGMs [16] has generally assumed prior knowledge of the class labels ai and in ([T]), or other covariates 
[T7H20]. Then, learning the parameters of ([T]) practically reduces to a simple logistic regression. However, the learning task is 
considerably more complicated if the latent class labels ai and are unknown and need to be inferred as hidden variables. On 
the other hand, a growing body of work is dedicated to the development of efficient algorithms for learning general stochastic block 
models [2TH25] including the hidden assignment of nodes into classes, but without the incorporation of node specific effects, i.e. a 
model specified by 

^(A. = 1\0) _ B..., 
P( = 0|^) 1 - B.^^^ ' 

This model is also referred to, with slight variations, as infinite relational model |26] or mixed membership stochastic block model [27] . 
The common thread among these is the goal to model network structure entirely in terms of groups of nodes. Some attempts to 
include the estimation of node specific effects have resulted in biased models [28] [29]. Within the framework of ERGMs, node and 
group specific properties have been combined in so called latent space models [30" "3T] where nodes are assigned a position in an 
abstract space and links form as a function of their distance. Such models are well motivated for social networks, where homophily 
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is a central mechanism of link formation and proximity in the latent space may be interpreted as similarity. Yet they are less general 
than stochastic block models being caught in the predicament of placing groups of nodes with similar interaction partners in close 
proximity while at the same time having to place them apart if these nodes are not densely connected among each other. Our 
approach facilitates parameter estimates and latent class inference in model ([T]) which combines node specific effects with the more 
general stochastic block models for group structure. 



B. Model Inference 

To describe the algorithm for estimating the parameters 0, we first write the likelihood of the entire observed network adjacency 
matrix A in terms of our model ([T]): 

c{e) = v{A\e) = Ylr{Di^ = A,^\e) (3) 

For a dyadic model, the likelihood factorizes into terms that involve parameters associated with only two nodes. 

Commonly used methods to estimate the parameters and hidden variables in such a model are to employ maximum likelihood 
(ML) techniques in the form of an expectation-maximization type algorithm or Monte Carlo sampling prefer a Bayesian 

approach, based on Maximum A Posteriori (MAP) estimates that does not incur the computational cost of Monte Carlo sampling 
while being less sensitive to initial conditions and more stable numerically than ML, especially as the parameters which maximize ^ 
may lie on the the borders of the admissible interval (0, 1). Furthermore, the MAP approach provides a natural Occam's razor as the 
posterior distributions of parameter estimates can only reduce in variance with the provision of more data, while the ML approach 
assumes point estimates or ^-functions for the posterior from the start. This is an important feature of the Bayesian approach as it 
provides a natural limit for the number of inferred classes and confidence levels in the assignments. Classes cannot be arbitrarily small 
if the posterior for the inter-class link preference B is to be localized. In contrast, under an ML approach the likelihood increases 
monotonically when more and hence smaller classes are used and model selection criteria, as in [20j, are needed. Finally, Bayesian 
techniques offer a principled way to incorporate prior domain knowledge for obtaining a more accurate approximate marginal posterior 
distribution V(Ok\A), where Ok represents one of the parameters ai^cn^ P^^r^ or Brs- 

A message passing or belief propagation algorithm provides a principled way to calculate approximate posterior marginal distributions 
[33]|34]. The starting point for this algorithm is a so-called factor- or dependency-graph, a graphical representation of the probabilistic 
dependencies between the variables (model parameters) we wish to infer from the data, and the individual factors that constitute the 
likelihood ([5]). Figure [l]^ shows this for the case of a bi-partite network, likelihood ^ and model ([T]). 

The algorithm proceeds by exchanging messages, conditional probabilities, between factors and variables connected in the depen- 
dency graph until convergence. Using the definitions: 

R^^.(Ok) = V(D^^ = A^^\Ok,A\A^^) and 

Q^^.(Ok) = r(ek\A\A^^), (4) 

one can interpret Ri^{Ok) {R-Message) as the likelihood of a single observed matrix entry given only the parameter Ok and 
all the data matrix except for entry Equally, Qi^{Ok) {Q-Message) is interpreted as the posterior probability distribution of 

parameter Ok given the entire data matrix except for entry Ai^. For the sake of notational economy, we have adopted to identify 
functions by their argument. It is to be understood that Ri^{ai) is a different function than Ri^{l3^) and not the same function 
Rii^{x) evaluated at the points at and /^^^ as should be clear from the definitions Q. 

Formally, we obtain the R-Message from Ai^ to Ok, by integrating out all parameters except Ok from a likelihood function 

R^^.{Ok) = ^V{D^^ = A^^\d,A\A^^)v{0\Ok\Ok,A\A^^)dd\Ok (5) 

Using the independence of given data entries Ai^ we can readily identify V(^Ai^\d^ A\Ai^^ with the of ([l]). Assuming the 

joint distribution V(0\A\Ai^^ factorizes with respect to every single Ok, one obtains the following closed set of equations: 

R^^.{Ok =x) = f P( = A^\0) n Q^^(Oe)dOe and 
^ i*k 

Q^^{Ok=x) oc V{Ok) n Rj^{Ok = x). (6) 

Although the factorization assumption may seem strong, it merely means that the Q-Messages V{Ok\A\Ai^) for any two variables 
Ok and Oi with k ^ I are assumed independent. Given that these distributions are conditioned on the entire data matrix except for 
one entry, the error we make using this assumption is considered negligible for large systems. The form of calculating Qi^{Ok = x) 
in ^ follows directly from Bayes' theorem and V{Ok) is the distribution we use to include prior information. These equations can 
be iterated until convergence after which we finally obtain the desired approximate marginal posterior distribution, for every single 
parameter, as: 

V{Ok\A) ocV{0k)Y\R^^{0k). (7) 
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FIG. 2: Attendance record of 18 women (rows) to 14 informal social events (columns) due to ethnographers Davis, Gardner 
and Gardner [35j. Black squares indicate attendance, a) Attendance matrix with posterior probability of class assignment for 
actors V{(j) and events V{t) as found by learning a standard stochastic block model Classification inferred divides events 
according to number of attendants and actors according to number of events participated in. The Inset shows the observed 
numbers of attendances do not agree well with the expectations due to model ([2]). b) The same attendance matrix as in a) 
but reordered due to the classification given in the original study indicated by dashed boxes 35 . Posterior probability of class 
assignments inferred using model ([T]) is almost perfectly compatible with the expert's classification. Including node specific 
popularity and activity parameters (3 and a allows to match observed numbers of attendances vs. expectations from model ([T]) 
as shown in inset. 



To illustrate these ideas, explicit update equations for the inference of the hidden class index at of node i appear below. Expressions 
for other parameters will be reported elswhere. With 

Qii^(ai)Qii^(Pfx)Qifx(Brs)daidf3^dBrs , (8) 
we can write for the R- and Q-Messages between Ai^ and ai\ 

Rifj.((Ti=r) = Y^XrsQif^ir^ = s) and 

s 

Q^^.(a^=r) oc r(a^=r)YlR^^.(a^=r). (9) 

The dependency graph greatly facilitates setting-up these update equations. Following the rules that R-Messages are always 
sent from factors to variables and Q-Messages from variables to factors; and that in R-Messages, we sum or integrate over the 
incoming Q-messages, while Q-Messages are proportional to the product of incoming R-Messages, we can write the equations based 
on the dependency graph. Figure [l^ shows a detail of[TJ\ focussing on factor to illustrate the messages involved in the 
calculation of Ri^(ai) sent to variable ai as in Generalizing the algorithm and update equations to undirected, uni-partite 

networks is straightforward. A more difficult task is the generalization to directed networks. Directed uni-partite networks are 
generally represented by an NxN adjacency matrix A. They may exhibit a feature that undirected networks cannot have: non-trivial 
reciprocity, i.e. the fact that the link from node i to node j is generally not independent of the link from node j to node i. Still, we 
can model such networks using dyads with i <j. But now, a pair of nodes can not only be connected or disconnected, but we 
have four different states a dyad can be in: it can be disconnected: Dij = (0,0), with a connection running only from node i to j: 
Dij = (1,0), or only from j to i: Dij = (0,1) or the connections between i and j may be reciprocated: Dij = (1,1). Hence, we 
have to incorporate a parameter that can model this reciprocity. Instead of one global parameter, here we allow reciprocities to differ 
depending on the latent classes: 



v{Dij = (1, i)\d) _ p,.,, r{D,, = i\§) v{D,r = \\e) 

(10) 



T{Dij = (0, 0)1^) (1 - ) V{Dij = Q\e) V{Dji = 0) 

p(A, = (i,o)|e) v{D,, = i\e) 
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Larger parameter values prs ^ (0, 1) increase the odds of finding reciprocated edges between nodes in classes r and s. Similar to Q, 
the parameters ai and pi model the tendency of node i to initiate and attract links, respectively. The model amounts to coupling 
two independent models of the type ([!]) via the reciprocity parameters. 

The likelihood of the observed data under such a model for given parameters is then given analogously to ^ by 

c{e)^v{A\e) = n^(A, = {A.j,A,,)\e) (ii) 

i<j 

and factorizes into terms involving only the parameters of two different nodes. Figure [Tfl shows the dependency graph for such a 
directed network. With the specification of the model and the dependency graph, the update equations are derived in the same way 
as above. 



II. RESULTS 



Next, we will demonstrate the impact of microscopic (node specific) effects on inferred mesoscopic latent class structure by 
comparing model ([T]) with the less expressive standard stochastic block model Q. For this, we use a dataset from sociology: the 
Southern Women. Then, we will show the importance of mesoscopic group effects to the interpretation of microscopic structural 
features by studying the motif distributions in the neural network of the nematode C. elegans. Last, we will determine the influence 
of both node specific and group specific effects on the accuracy of predicting new links in networks. To this end, we compare model 
([T]) with both a less expressive model and a more expressive model in terms of classification accuracy and predictive power on the 
network of gene-disease associations from the Online Mendelian Inheritance in Man (OMIM) database. 

Southern women: This classic bipartite data set is due to ethnographers Davis, Gardner and Gardner [35]. A 18 x 14 matrix records 
the attendance of 18 women in southern Alabama to 14 informal social events over the course of a nine month period in the 1930s. 
The authors' aim was to study how an individual's social class influences her pattern of informal social interaction. Based on intuition 
and experience in the field, but without formal analysis, the authors suggested the existence of two latent classes of 9 women each, 
with only little overlap in the attendance at events. Over the years, the data has become a standard test case of network analysis 
algorithms, a meta-analysis of which can be found in [36j. We are interested in whether an inference based approach can assert the 
presence of latent classes and whether the class assignments found correspond to those suggested by the field experts. 

First, we learn the standard stochastic block model with two latent classes for both actors and events, i.e. model and only 
estimate class membership ai, and preference matrix Brs- Figure[2^ shows the data, with rows and columns of the attendance 
matrix reordered such that events/actors predominantly assigned to the same class are adjacent. The resulting block model is in 
contrast to findings of the original authors [35j. Events seem divided according to number of participants, i.e. popularity, while actors 
seem divided according to number of events participated in, i.e. activity. The expert classification due to social class is not correctly 
captured when trying to model the network through group effects alone. The reason for this failure is that under model Q, the 
degree distribution for members of the same latent class is assumed to be Poissonian. The expected degree is the same for each 
member of a class. The inset in figure shows that this assumption cannot capture the observed degree distribution. Since the 
standard stochastic block model does not model node degree independently of class preference; variance in degree distributions of 
both actors and events confuses the inference of group membership. 

In contrast, the inset in figure |2Jd shows the expected degree vs. the observed degree when including activity and popularity 
parameters in the model as in ([T]) and allowing for two classes. Now, the observed degree distribution can be accounted for. The 
introduction of activity and popularity parameters has also dramatic effects on the latent classes inferred. Figure |2Jd again shows 
the attendance matrix, but this time, rows and columns are ordered as given in [35j and the authors' assignment to social class is 
indicated by dashed boxes. The experts' classification matches almost perfectly that inferred using model ([T]). We can see that events 
such as 8 and 9 which are attended by most actors receive high (3 values and thus have very little discriminative power. Also, actors 
who are very active and occasionally participate in events predominantly frequented by actors from the other group such as Mrs. N. 
F. can still be assigned with high probability to a class, despite conflicting evidence in their participation record. Using model ([T]) 
effectively allows one to decouple the preference effects of a group of actors for a group of events from global effects that contribute 
to the variance in node connectivity. 

Caenorhabditis elegans: In our second application, we explore the extent to which a dyadic model may explain the distribution of 
small sub-graphs, termed motifs, in a network. Motifs have received considerable attention as possible entities of network formation, 
i.e. building blocks larger than single edges. Their distribution relative to random null models has been suggested to characterize 
entire classes of networks [8J. The over/under-representation of certain motifs with respect to random null models is often attributed 
to possible evolutionary pressures due to a motif's potential influence on the performance of the network's function [38 39j. 

We study the distribution of all 16 possible 3-node motifs in the 279 neuron chemical synapse network of C. elegans [37j. The null 
model commonly used to assess whether a particular motif is under- or over-represented in a network is generated by randomizing the 
original network conserving only microscopic structural features, i.e. the number of incoming, outgoing and reciprocated links at each 
node is preserved. All other structural features and correlations are removed by the randomization. Figure [3jD shows box-plots for 
motif counts in 1000 such random networks and in comparison the actual count of the 16 motifs in the chemical synapse network of 
C. elegans normalized to the mean count found in the set of null models. We can see that using this null model, 11 of the 16 motifs 
are strongly over/under-represented and hence would qualify as possible starting points for further research on putative functional 
relevance. 

However, the standard null model also removes all mesoscopic structures, in particular structure due to groups of more than 
three nodes. A dyadic model such as (10) lacks any parameter for three-node motifs but can generate an ensemble of null models 
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FIG. 3: Motif counts in the synapse network of C. elegans compared to two random null models, a) Adjacency matrix of the 
observed neural network Adjacency matrix of a typical realization of a link randomized version of the original data 

and resulting Z-score statistics of motif counts. Counts in the original data (red x) are compared to box plots of counts in 1000 
link randomized null models. Strong deviations are found at 11 of the 16 motifs. Since the link randomized null models retain 
only node specific features, i.e. the numbers of incoming, outgoing and reciprocated links at each node, the cannot capture the 
apparent mesoscopic structure in the original network and henc e m ay over-estimate the statistical significance of some motifs, 
c) Adjacency matrix of a typical network generated from model ( 10 ) with both node specific as well as class specific parameters 
estimated from the original network. 15 classes were used in this example. Using 1000 networks generated from model (10) as 
a reference ensemble, the Z-score statistics show mild deviations only at 3 of the 16 motifs. This indicates that class structure 
may offer a more parsimonious explanation for the observed motif distribution. 
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FIG. 4: Classification accuracy and predictive power of network models ([T|, Q and that by Newman/Leicht (NL) [23 . a) 
Overlap of an expert classification of diseases in the Diseasosome-Network 40 and that inferred using models and the data of 
known gene-disease associations recorded in the Online Mendelian Inheritance in Man (OMIM) database by Dec. 2005. Measure 
of overlap is normalized mutual information (NMI) [41 . b) Prediction accuracy at 16 classes for confirmed associations added 
to the OMIM database between Dec. 2005 and Jun. 2010. For each model, a candidate list of associations is obtained by 
sorting all possible associations in descending order according to their probability under that model with parameters estimated 
from the Dec. 2005 data. We plot which fraction of actually confirmed associations is found in the corresponding top fraction 
of the candidate list. Entries due to new variants of a previously recorded association are listed as "repeated associations" while 
genuine new associations are reported as "new associations". For example: In the top 1% of any candidate list, we expect to 
find 1% of new associations due to chance alone. We do find 15% of all confirmed new associations if the list was due to model 
([2|, 20% if the list was due to the NL model and 30% if the list due to model ([T]). See text for details. 



that nnatches the observed network in ternns of the observed node specific degrees as well as with respect to nnesoscopic structural 
features. Such nnesoscopic structure inevitably exists as neurons are located in different sonnatic regions and synaptic connections 
between closely located neurons are nnore likely than between distant ones [42J. Neurons are also aggregated in different ganglia 
nnaking intra-ganglia connections nnore likely than inter-ganglia synapses. Furthernnore, they serve different functions that influence 
their connectivity. For exannple, stinnuli nnay be processed in a sensory neuron - interneuron - nnotor neuron cascade. The latent 
classes we infer fronn the data using (10) can be explained using a connbination of these factors (see suppl. nnaterial). For instance, 
Sonne classes connbine nnotor neurons only, sonne connbine both sensory, inter- and nnotor-neurons which are all located in a particular 
sonnatic region, while others can be explained post-hoc as lying in the sanne ganglia. More innportant than the interpretation of these 
classes is whether a dyadic nnodel that assunnes all pairs of nodes as conditionally independent can account for the observed three 
node nnotif-counts in the network. 

Figurejs]: shows the box-plots of nnotif counts in 1000 networks generated fronn ( 10 ) allowing for 15 different classes of neurons and 
using the paranneters estinnated fronn the original network, again nornnalized to the nnean count. The connparison with the nnotif-count 
in the C. elegans network now shows that only 3 out of 16 nnotifs cannot be explained by the null nnodel. This result is rennarkable 
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as it underscores the importance of group specific effects in modeling complex networks. The fact that a simple dyadic model can 
explain a large portion of the three-node statistics in the observed data is a strong corroboration for our claim that latent classes of 
nodes are important determinants of network structure. Furthermore, it offers a very parsimonious explanation of motif statistics in 
this network and a more conservative estimation of their statistical significance. 

OMIM: As a last example, we study the influence of classification accuracy on the predictive power of probabilistic models using 
a bi-partite network known as the human "Diseasosome-Network" [40j. It represents known associations between genes and diseases 
recorded in the OMIM database [43]. The network was first published in 2005 and we focus on the analysis of the largest connected 
component involving 516 different diseases and 903 different genes connected by 1550 different associations known in 2005 [40j. The 
original publication provided an expert classification of the diseases into 22 types. The type of disease is predominantly based on the 
tissues and organs involved (such as bone, connective tissue, muscular, dermatological, hematological, renal, etc.) or based on the 
affected system (such as skeletal, cardiovascular, immonological, metabolic or endochrinal, etc.) 

To what extent does such a classification overlap with one inferred from a network of common genetic causes? We compare model 
([T]) with the less expressive standard stochastic block model ([T]) and a more expressive model due to Newman and Leicht (NL) ^29] . 
The latter includes both individual and group effects as in ([T]), but instead of a single parameter for the overall activity or popularity 
of a node, it features one such parameter per latent class, i.e. it models activity or popularity with respect to each class of nodes. 

We compare the overlap between the expert classification of diseases and the one found algorithmically, based on the gene-disease 
association network alone. We restricted ourselves to using the same number of classes for both genes and diseases. The comparison 
of models ([T]), NL and the standard stochastic block model ^ is shown in figure [4^. As expected from the earlier discussion, 
neglecting individual node effects as in model ^ reduces the overlap with an expert classification compared to model Q. But, 
interestingly, the same applies when including gene-specific effects for every class of diseases and disease-specific effects for every 
class of genes as in the NL model. Too many explanatory variables per individual node seem to reduce the detection quality of latent 
classes. 

Since 2005, the OMIM database has been steadily growing and 292 new associations between those 516 genes and 903 diseases 
had been added until June 2010. Using the data from 2005 as a training set and these new additions as a test set, we compare the 
predictive power of the different models for future associations. New entries to OMIM comprise both new variants of already known 
gene-disease associations (repeated associations) as well as genuine new associations of genes with diseases that were not linked 
previously. Hence, the data offers the opportunity to differentiate predictive power with respect to these two types of entries. Using 
the parameters estimated from the 2005 data set for each model Q, NL and Q, we calculate the probability for association of each 
gene i with each disease fi as V{Di^\0). Then we sort these probabilities in descending order and hence obtain a candidate list for 
new or repeated associations. For instance, in the case of models with 16 classes, figure [4^ now shows how far we have to go down 
the candidate list to find a certain fraction of the associations that were added to the database over the course of 4 1/2 years. 

Variants of already known associations seem to be added approximately randomly to the database as models Q, NL and ^ 
all perform close the random expectation for these repeated associations. For the genuinely new associations, however, we observe 
that all models strongly deviate from the random expectations. In particular ([T]) outperforms both NL and Q, with the latter two 
performing similarly. 

Comparing figures[4^ and[4jD, we conclude that the standard stochastic block model may be too simple to capture the biologically 
relevant network structure as seen from a low overlap with an expert disease classification and low predictive power. With the inclusion 
of node specific effects, the NL model is more flexible in capturing the observed network leading to a higher overlap with the expert 
classification of diseases. It does, however, not generalize well as its predictive ability is still lower than that of the ERGM model 
([T]). The latter appears to provide the best compromise between flexibility and parsimony as it combines excellent ability to represent 
biologically relevant structures with high classification accuracy and a parsimonious inclusion of node-specific effects, leading to the 
best predictions. 



III. DISCUSSION 

We have presented an efficient, distributive algorithm that successfully estimates the parameters and latent group assignments of 
an exponential random graph model including both node specific and group specific properties. We have shown that including node 
specific effects in the estimation of latent classes leads to improved recovery of class assignments by domain experts. Additionally, 
we have shown that using a simple dyadic model, a large part of the triad statistics in networks may be explained, shedding new 
light on the discussion of motif distributions in complex networks. We expect our results to stimulate a discussion on the use of 
appropriate null models in the analysis of sub-graph distributions and their universality for certain classes of networks. Finally, we 
have explored the predictive power of the model to identify new gene-disease associations, using the OMIM database. Through these 
specific examples, we have demonstrated that node specific and group specific properties should be both incorporated when inferring 
and modeling structural features in complex networks. 
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